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Abstract 

We present the characterization and calibration of the slow-scan observation mode 
of the Infrared Camera (IRC) on-board AKARI. The IRC slow-scan observations were 
operated at the S9W (9 jttm) and L18W (18 jttm) bands. We have developed a toolkit 
for data reduction of the IRC slow-scan observations. We introduced a "self-pointing 
reconstruction" method to improve the positional accuracy to as good as 1". The 
sizes of the point spread functions were derived to be ~ 6" at the S9W band and 
~ 7" at the L18W bands in full width at half maximum. The flux calibrations were 
achieved with the observations of 3 and 4 infrared standard stars at the S9W and 
L18W bands, respectively. The flux uncertainties are estimated to be better than 
20 % from comparisons with the AKARI IRC PSC and the WISE preliminary catalog. 
Key words: infrared: general — techniques: image processing — space vehicles 



1. Introduction 

AKARI is the first Japanese infrared astronomical satellite (Murakami et al. 2007), which 
was launched on 2006 February 21. AKARI was brought into a sun- synchronous polar orbit at 
an altitude of 700 km. AKARI has a 68.5 cm diameter cooled telescope with the two scientific 
instruments of the Infrared Camera (IRC; Onaka et al. 2007) for 1.8-26 fim and the Far-Infrared 
Surveyor (FIS; Kawada et al. 2007) for 50-180 /im. One of the major observation programs of 
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AKARI is an all-sky survey at the mid- to far-infrared wavelengths with 6 photometric bands to 
provide second generation infrared catalogs. In addition, AKARI has the capability for imaging 
and spectroscopy in the pointed observation mode. 

The IRC is mainly designed for deep imaging and spectroscopy over 1.8-26 /im with 
three independent channels: NIR (1.8-5.5 /im), MIR-S (6-12 /im), and MIR-L (12-26 /im). All 
the channels have filter wheels, which hold three filters and two spectroscopic dispersers. The 
field-of-view (FoV) is 10' x 10'. The MIR-S and MIR-L channels have infrared sensor arrays 
of 256 x 256 pixels. The pixel sizes for MIR-S and MIR-L are 2"34 x 2"34 and 2"51 x 2"39, 
respectively. The FoVs of MIR-S and MIR-L are separated by 20' perpendicularly to the 
AKARI's scan direction. 

Although IRC was originally designed for imaging and spectroscopy in the pointed ob- 
servation mode, Ishihara et al. (2006) developed an additional function, "scan mode", to carry 
out the All-Sky Survey with the MIR-S and MIR-L channels. In the "scan mode" , we used 2 
out of 256 rows in the sensor array for continuous and non-destructive readout. The data were 
sampled at every 0.044 sec, and the pixels were reset at every 13.5 sec (or every 306 samplings) 
to discharge the photo-current. The output signals of every four adjacent pixels along each 
row were binned on board to reduce the data amount, so the effective pixel size in the cross- 
scan direction was about 10". The details of the IRC "scan mode" and the all-sky survey are 
described in Ishihara et al. (2006) and Ishihara et al. (2010). 

2. IRC Slow-scan Observations 

The IRC slow-scan observation mode was designed for large-area (up to 10' x l.°5 in one 
pointed observation) mapping observations with sufficient sensitivities, using the "scan mode" . 
The IRC slow-scan observations are provided in the Astronomical Observation Templates 
(AOTs) of IRC11 and IRC51. In the slow-scan observations, the telescope scanned along 
one or two round trip paths around the target object to take a redundant dataset. More details 
of the telescope operations in the slow-scan observations are described in Kawada et al. (2007) 1 . 

The basic properties of the slow-scan mode are the same as in the All-Sky Survey: (1) 
the data of the S9W (9 /im) and L18W (18 /im) bands were obtained with the MIR-S and 
MIR-L channels, respectively 2 , (2) the data of only 2 rows in the 256 x 256 pixel array in each 
channel were acquired, and (3) the data were sampled at every 0.044 sec with non-destructive 
readout. On the other hand, there are three points different from those in the All-Sky Survey 
as follows: 

1. The pixels were reset at every 306 and 51 samplings to discharge the photo-current for 

1 IRC and FIS were simultaneously operated in the slow-scan mode. 

2 Only one set of observations was carried out with the Sll (11 /im) band with the MIR-S channel (Ishihara 
et al. 2007). 
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IRC11 and IRC51, respectively. 

2. The telescope scanned the sky with a much slower speed (8, 15 or 30" sec" 1 ) than that of 
the All-Sky Survey (216" sec" 1 ). 

3. The output signals of every four adjacent pixels along each row were binned on board, 
so that the effective pixel size in the cross-scan direction is 10" for IRC11 (same as the 
All-Sky Survey), while no binning was applied for IRC51 so that the pixel size is 2"5. 

Since the array operations of IRC11 are identical with the All-Sky Survey and most of 
the slow-scan observations were performed with IRC51, we describe IRC51 observations in this 
paper. 

3. Data Packages of the IRC Slow-scan Mode 

Data are stored in dedicated format FITS file, called Time-Series Data (TSD), which 
was originally developed for the FIS (see Verdugo et al. 2007), for the data package of the IRC 
slow-scan observations to handle the data easily. The TSD is a binary FITS table and consists 
of a header part and arrays of data records. One record consists of the IRC instrument raw data 
and necessary information from other house keeping (HK) instruments, as well as positional 
information given by the ground attitude determination processing. Since the sampling rates of 
these data are different from each other, the data should be interpolated to synchronize exactly 
with the clock of the detectors. Two data sets for the MIR-S and MIR-L channels of the IRC 
are created from one pointed observation. 

4. Data Reduction Processes 

We have developed a dedicated software for reduction of the IRC slow-scan data, ARIS 
(AKARI data Reduction tools for the IRC Slow-scan). ARIS is written in Interactive Data 
Language (IDL). Most of the ARIS processes handle the TSD format data. We describe indi- 
vidual data reduction processes in the following sub-sections. 

4-1. Basic calibrations 

In the first step of data processing, we applied the following basic calibrations. First, 
we corrected the anomalous behaviour of the detector output, which is seen after the reset. 
Second, we corrected the non-linearity between incoming photons and output signals. Third, 
we differentiated the data, because the output signals from the detectors are time integral 
values. Fourth, we subtracted the dark signals. Then, we applied the flat fielding. Finally, we 
masked bad data, such as pixels which were masked out for the slit spectroscopy, pixels just 
after the resets, and saturated pixels. 
Reset anomaly correction 

There is an anomalous behaviour of the detector output which persists for a few seconds 
after every reset (the reset anomaly). Figure 1 shows an example of the anomaly. Since the 
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Fig. 1. Example of the reset anomaly correction. The plus and cross symbols are the output signals of a 
pixel from observed sky and corrected data, respectively. 

Table 1. Parameters of the reset anomaly correction. 





MIR-S 


MIR-L 


ao 


3.098 


2.721 


Oi 


-0.1045 


-0.08774 


«2 


72.31 


4.973 


«3 


-0.4321 


-0.1087 



offset level of the detector output is sensitive to the temperature (Ishihara et al. 2003), the 
anomaly can be well described by the temperature drift caused by the tiny reset current to 
discharge the stacked photo-electrons. The reset anomaly is corrected as 

sP(t) = S?\t) + R{t), (1) 

where S^, S^, and R are the raw data, the corrected data, and the correction offset, respec- 
tively. The suffix i indicates the pixel number and t is the sampling number from the last reset. 
The offset is given by 

R(t) = a x exp(a! x t) + a 2 x exp(a 3 x t 2 ). (2) 

The parameters we adopted are listed in Table 1. 
4-1.2. Non-linearity correction 

The IRC detectors have non-linearity due to the decrease of the bias voltage during 
integration. This non-linearity was measured in the laboratory prior to the launch. Figure 2 
shows examples of raw signal output vs expected signal, where the expected signal was derived 
by fitting a linear function to the raw data of ADU < 15000 with a linear function. We adopt 
the correction factors in a polynomial form as 
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Fig. 2. Plots of raw signal versus expected signal for the MIR-S (left) and MIR-L (right) channels. The 
green curves are the best-fit polynomial functions given by equation (3). The broken lines indicate the 
best-fit linear function for ADU<15000. 
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Fig. 3. Deviations from the best-fit polynomial functions (green dots) for the MIR-S (left) and MIR-L 
(right) channels. The red dots are the deviations from the best-fit linear functions for ADU<15000. 



sP^J2L x x(S^r, 



(3) 



x=l 



where L x is the coefficient for the term of degree x. We adopt n = 7 for both the MIR-S and MIR- 
L channels, and the best-fit parameters are listed in Table 2. Figure 3 shows the deviations from 
the best-fit polynomial functions, indicating the polynomial fitting is as accurate as 5 % for the 
output signals less than 30000 and 35000 ADU for the MIR-S and MIR-L channels, respectively. 
Figure 3 also shows that the physical detector saturation occurs around 40000 ADU for both 
the channels. 
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Table 2. Best-fit parameters of the non-linearity correction. 





MIR-S 


MIR-L 


Li 


0.9977 


0.9996 


L-2 


2.559 x 1(T 6 


-2.906 x 10~ 7 


L, 


-4.029 x 10- 10 


5.278 x 10~ 10 


Li 


-7.599 x 10~ 15 


-1.132 x 10~ 13 


L 5 


3.481 x 10~ 18 


8.722 x 10~ 18 


L 6 


-1.387 x 10~ 22 


-2.659 x 10~ 22 


L 7 


1.722 x 10~ 27 


2.913 x 10~ 27 



4-1.3. Differentiation 

Since the raw signal from the detector in scan observations are time-integrated data, we 
differentiate it with respect to sampling number t as 

S?\t) = S?\t)-Si 2 \t-l). (4) 

4-1. 4- Dark subtraction 

The slow-scan observation mode took dark frames during the manoeuvre at the beginning 
and ending of each scan observation. We calculate the dark current per pixel by taking the 
median value of all dark frame data, and subtract it from the pixel value, assuming that the 
dark level was almost constant in one scan, 

Sf^-A, (5) 

where D is the dark signal. 
4-1.5. Flat fielding 

We construct the flat-field data from ~200 observations of the south ecliptic pole (SEP) 
region by taking the median value for each pixel. The correction is applied as 

Sf = S\ 4) /F u (6) 

where F is the normalized flat correction factor (see Figure 4). 
4-1.6. Masking bad pixels 

We set flags on invalid data such as 

1. pixels which were masked out for the slit spectroscopy, 

2. data of the first 2 samples after each reset because of significant reset anomaly, and 

3. saturated data which have ADU>40000 in the raw signal output. 

These data are not used to create images. 
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Fig. 4. Flat data for the S9W (bottom) and L18W (top) bands normalized by the mean values of pixels 
25-256 in row#l and 281-512 in row#2 for the S9W band and 2-236 in row#l and 258-492 in row#2 
for the L18W band. 

4-2. Image creation 

After the basic calibrations, we proceed to image creation. In this process, we use a 
GCF (Gridding Convolution Function; Sawada et al. 2008) to re-grid the sampling data onto 
a regular grid, since the data acquisition was spatially irregular in the sky plane. The default 
GCF is a Gaussian function with full width at half maximum (FWHM) of 1 pixel, where the 
pixel sizes are 2."34 and 2."51 at the S9W and L18W bands, respectively. To abandon charged- 
particle-hit data, we perform 5-cr clipping within 1 pixel radius 3 . We note that threshold lower 
than 3-a will diminishes the peak intensity of a true source in some degree. 

4-3. Self-pointing reconstruction 

In making images, we need positional information of every pixel at every sampling. There 
are two kinds of information in TSD; AOCU (the Attitude and Orbit Control Unit) and G-ADS 
(the Ground-based Attitude Determination System) (see Verdugo et al. 2007). However, their 
positional accuracy is typically 10" in peak-to-peak, (see the left image in Figure 5). Although 
these accuracies satisfy the requirement for the attitude controle system of AKARI of 30" 
absolute error, it is not accurate enough compared to the pixel size of 2"5 and 5" in IRC51 
and IRC11 observations, respectively. Actually, there appear many fake double stars owing to 
the positional differences between forward and backward scans in each round trip. Therefore, 
ARIS has a function to make time-dependent correction of scanned position based on the TSD 
by comparing the positions of detected point sources with those in the reference catalogs. This 
process is called self-pointing reconstruction (self-PR). 

The self-PR procedure is performed as follows. 

1. We make two or four initial images of each individual one way scan in each round trip 



3 Usually there are more than 20 data in 1 pixel radius, although the number varies with the scan speed. 
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using either AOCU or G-ADS position data (see the left panel in Figure 5). Usually, we 
use the S9W data because we need as many point-like sources as possible in the images. 

2. We extract point sources and store their observed time and positions. To avoid extracting 
fake sources, we use a high threshold level of 5-er. 

3. We do cross-identification between the detected sources in our images and the point sources 
listed in the reference catalogs. As a reference catalog, we usually use the Two-Micron All 
Sky Survey (2MASS) Point Source Catalog, whose positional accuracy is as good as 1". 
Since the 2MASS sensitivity is as good as 14 mag in the K$ band, most of the AKARI 
sources are easily identified in the catalog. 

4. We evaluate the differences between our positions and those in the catalog. Here, the 
positional differences are measured in the in-scan and cross-scan directions. 

5. We take mean values of the differences in a given time interval (30 sec as a default value), 
and then, connect the values by line segments to get the overall trend (see the middle 
panel in Figure 5). The computed trends are output as a text file as time vs. positional 
difference format. 

6. We revise the position information in TSD with the computed trends and make the final 
image from all forward/backward scans. 

The positional uncertainty in the final image is as small as 1" (see the right image in 
Figure 5). In addition, we can also correct the positions for the L18W and the FIS slow- 
scan data, which were observed simultaneously, using the trend file made from the S9W data. 
The application to the FIS data was demonstrated in Ikeda et al. (2012). We note that the 
performance of the self-PR method becomes maximum when sufficient number of bright S9W 
sources are distributed over the scanned area. 

5. Flux Calibration for Point Sources 

5.1. Observed standard stars and data processing 

For the absolute flux calibration, we observed 4 and 3 infrared standard stars (Cohen 
et al. 1999) for the S9W and L18W bands, respectively (Table 3). These observations were 
carried out in the AOT IRC51. 

The data were processed by ARIS, and the final images were created using the self-PR 
method. We applied aperture photometry to the final images. We used the DAOPHOT-Type 
Photometry Procedures of GCNTRD and APER in the IDL Astronomy User's Library (version 
2007 May; Landsman 1993). The radius of the aperture was 7.5 pixels and the sky level was 
estimated from the annulus between 7.5 and 12.5 pixel radii. 

5.2. Estimation of the in-band flux density 

The in-band flux density of each band at the effective wavelength, /^ uotcd (Aj) is calculated 
by the following equation: 
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Fig. 5. Example of the self-pointing reconstruction, (left) The 9 fxm image from the first forward-scan 
using the satellite telemetry and a magnified image of a point source. The pixel size is 2."34. The crosses 
are the positions of 2MASS point sources, (middle): Time variations of the positional differences along 
the cross-scan (top) and in-scan (bottom) directions, respectively. The thick curves indicate the overall 
trends of the time variations of the positional differences, (right): The final image based on the data 
positionally corrected by the self-PR method. 



Table 3. Observed standard stars. 



Expected flux density [Jy]t Observation 



Star 


Sp. type 


S9W 


L18W 




date (UT) 


ID 


par am* 


BD+62 1588 


K5 III 


4.9966X10" 1 


1.1341x10- 


i 


2007-03-13 23:59:21 


5124045 1 


i;N;15 












2007-03-14 01:38:43 


5124046 1 


i;N;30 












2007-03-16 00:00:36 


5124047 1 


i;L;15 












2007-03-16 01:39:57 


5124048 1 


i;L;30 


HD 42525 


AO V 


2.7792X10" 1 


5.8763x10" 


2 


2007-05-02 09:00:33 


5124066 1 


i;N;30 












2007-05-03 01:34:32 


5124065 1 


i;N;15 












2007-05-04 23:57:43 


5124067 1 


i;L;15 












2007-05-05 01:37:07 


5124068 1 


i;L;30 


HD 46819 


K0 III 


4.4501 xlO" 1 


9.6344x10" 


2 


2007-06-06 00:46:43 


5124108 1 


a;N;15 












2007-06-06 02:26:08 


5124109 1 


a;N;30 












2007-06-17 01:56:29 


5124110 1 


i;L;15 












2007-06-20 02:51:59 


5124111 1 


i;L;30 


BD+66 1060 


K2 III 


1.1766X10" 1 


2.6382x10" 


2 


2007-08-02 04:29:25 


5124130 1 


a;N;30 



f Flux density for each band estimated from equation (7). 

\ The first characters 'a' and T indicate one and two round trips, respectively. The second parameter shows 
that the target source is observed with the MIR-S ('N') or MIR-L ('L') channel. The last numeral is the scan 
speed in arcsec sec -1 . 
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Fig. 6. Estimated flux density versus observed ADU of the standard stars for the S9W (left) and L18W 
(right) bands. A broken line in each plot represents the least-square fit to the data. 



Table 4. 


Conversion factor 


cf [Jy ADU- 1 ] 


Band 


cf 


Error 


S9W 


2.518 x 10~ 4 


3.228 x 1(T 6 


L18W 


5.291 x 10" 4 


1.365 x 1CT 6 



(.quoted 



f£Ri(\)\f x (\)d\ 
/ A A ! ;(A,,/A) J R J (A)ArfA' 



(7) 



where f\ (A) is the Spectral Energy Distribution of a standard star (provided by M. Cohen) 
and Ri(X) is the spectral response function (including the transmission of the optics and the 
response of the detector) of the band i. The in-band flux density becomes equal to a true value, 
only if f\ oc A -1 . This is the convention adopted by IRAS, COBE, ISO, Spitzer/IRAC, and the 
other AKARI observation modes (Tanabe et al. 2008; Ishihara et al. 2010). 

5.3. Conversion factor 

The relations between the intensities in ADU from the images and the calculated in-band 
flux densities in Jy of the standard stars are shown in Figure 6. Straight lines were fitted to 
these data and the slopes of the fitted lines provide the conversion factors cf in Jy ADU -1 , as 
listed in Table 4. 

We note that the calibrations for diffuse emission have not been established yet. 
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S9W L18W 



Fig. 7. PSF images of the S9W (left) and LJSW (rwjffa) bands. The pixel scale is 2!'34 and 2"51 for 
the S9W and L18W, respectively. The black and white contours in each panel represent the 10 % and 
50 % levels of the peak value, respectively. 

6. Performance of IRC Slow-scan Observations 

We evaluate the performances of the IRC slow-scan observations by detected point 
sources, such as detection limits, positional accuracies, and flux uncertainties, on the basis 
of the 146 pointed observations toward the Chamaeleon region (>30 deg 2 ). 

At first, we constructed images of the S9W and L18W bands using ARIS, applying the 
self-PR method. We have succeeded in extracting 2974 and 492 point-like sources from the 
S9W and L18W images, respectively, and 246 sources were detected in both the bands. We 
created point spread functions (PSFs) with 25 brightest and relatively isolated stars, excluding 
saturated ones, in the S9W and L18W images toward the Chamaeleon I molecular cloud using 
the 'GETPSF' procedure in the IDL DAOPHOTO package. Figure 7 shows the derived PSFs; 
the FWHM sizes of the PSFs are ~ 6" and ~ 7" for the S9W and L18W bands, respectively. 

6.1. Sensitivity 

Figure 8 shows histograms of the extracted sources by the Chamaeleon observations. 
The lower limit in S9W and the turnover in L18W are thought to correspond to the detection 
limits of the AKARI IRC slow-scan observations; 9.5 mag for the S9W band and 6.5 mag for 
L18W. These limits are superior to those of the AKARI All-Sky Survey of 0.05 Jy (7.6 mag) 
and 0.09 Jy (5.3 mag) at the S9W and L18W bands, respectively (Ishihara et al. 2010). Note 
that the sources brighter than 2 mag at the S9W band are thought to be saturated. 

6.2. Positional accuracy 

We checked the positional accuracy of the detected sources by comparing with the 
2MASS Point Source Catalog. We searched 2MASS counterparts with positional difference 
of 3" (see Ita et al. 2010) and found that 3074 out of 3220 AKARI Chamaeleon sources have 
counterparts. Figure 9 shows the positional differences between our sources and the counter- 
parts. The mean and standard deviation values of the differences are calculated as 0."87±0."58, 
which is negligible for the full IRC resolution of 2''5. We note that we cannot see any corre- 
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Fig. 8. Histograms of the magnitudes of the extracted sources in the S9W {left) and L18W {right) 
bands. 
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Fig. 9. Histograms of the angular distances between the AKARI S9W and 2MASS sources. 

lations between positional differences and flux densities. There are 146 sources (126 for S9W 
and 22 for L18W) without 2MASS counterparts. Since about 90 S9W sources without 2MASS 
counterparts are as faint as >8.5 mag, we cannot exclude the possibility that these faint sources 
are fake ones considering the detection limit of 9.5 mag at the S9W band. Furthermore, since 
most of other sources are located near the edge of the images, they are also thought to be 
fake ones. On the other hand, the absence of 2MASS counterparts for 22 L18W sources were 
thought to be caused by the limitation of the self-PR method; since the self-PR method was 
applied using bright S9W sources, we could not improve the positional accuracy for regions 
where no bright S9W sources were detected (see §4.3). Indeed, 16 out of 22 sources have 
2MASS counterparts within a positional difference of 5". 
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Fig. 10. Flux density ratios of the slow-scan to the IRC PSC as a function of the slow-scan flux densities 
for the S9W {left) and L18W {right) bands. 

Table 5. Mean values (MEAN) and standard deviations (STDDEV) of the flux ratio of the slow-scan to the All-Sky Survey. 



Band 


MEAN 


STDDEV 


n star 


S9W 


0.945 


0.178 


613 


(>0.1[Jy]) 


0.986 


0.160 


423 


L18W 


1.027 


0.195 


150 



-k 'n_star' indicates the number of the stars used in the estimation. 
6.3. Flux uncertainty 

We compared the flux densities of the detected sources with those in the AKARI IRC 
PSC to verify the flux calibration of the IRC slow-scan observations. There were 613 and 150 
sources detected at both the slow-scan and All-Sky Survey in the S9W and L18W bands, 
respectively. Figure 10 shows the ratio of the slow-scan to PSC flux densities as a function of 
the slow-scan flux density. It seems that the ratios for the faint (<0.1 Jy) S9W sources drops 
systematically down to ~0.5-0.7. This disagreement is thought to be caused by the 'flux-boost' 
effect in the IRC PSC. It is explained as follows; if a source is as faint as the noise level, we can 
detect it only when the signal is unevenly bright. Actually, the completeness of the IRC PSC 
is only ~50 % at /(9/xm) = 0.1 [Jy] (see Figure 20 and Table 3 of Ishihara et al. 2010). The 
mean values and standard deviations of the ratios are listed in Table 5. 

We also searched for the counterparts in the first preliminary release of the WISE (the 
Wide-field Infrared Survey Explorer; Wright et al. 2010) catalog to verify AKARI flux densities. 
WISE surveyed the entire sky at 3.4, 4.6, 12 and 22 /im bands (which are referred to as Wl, 
W2, W3, and W4) with spatial resolutions of 6" at W1-W3 and 12" at W4. The WISE 
preliminary catalog covers our surveyed region. We found that 3088 out of 3220 sources have 
WISE counterparts within 6" searching radii. Here, we considered only the WISE sources with 
SNR > 10 in the 12 /im bands. Figure 11 shows the S9W-W3 vs. S9W and the L18W-W4 
vs. L18W color-magnitude diagrams. In the diagrams, excess sources such as young stellar 
objects, asymptotic giant branch stars, and galaxies are located on the right-hand side. For most 
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Fig. 11. Color-magnitude diagrams of S9W-W3 (WISE 12 /mi) vs. S9W [left) and L18W-W4 (WISE 
22 /tm) vs. L18W (right). 

sources, AKARI and WISE measurements are in good agreement to each other considering the 
differences in the band passes; the mean colors are 0.12 and 0.037 at S9W- W3 and L18W- W4, 
respectively. 

7. Summary 

We have developed a toolkit, ARIS (AKARI data Reduction tools for IRC Slow-scan), 
and described in detail the data reduction processes. Using the self-PR method in ARIS, the 
positional accuracy of final images reaches as good as 1", smaller than the pixel size of IRC of 
2"5. The PSF sizes are ~ 6" and ~ 7" in FWHM for the S9W and L18W bands, respectively. 
We also performed flux calibrations using infrared standard stars. We further checked the 
uncertainties of flux densities comparing with the AKARI IRC PSC and the WISE preliminary 
catalog, and found that the uncertainties are less than 20 %. 

This work is based on observations with AKARI, a JAXA project with the partic- 
ipants of ESA. We gratefully acknowledge F. Usui, C. Yamauchi and all the members of 
the AKARI project for their support to this study. This study uses data products from the 
Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and 
the Infrared Processing and Analysis Center /California Institute of Technology, funded by 
the National Aeronautics and Space Administration and the National Science Foundation. 
WISE is a project of University of California, Los Angeles, and Jet Propulsion Laboratory 
(JPL) /California Institute of Technology (Caltech), funded by the National Aeronautics and 
Space Administration (NASA). This research has made use of the SIMBAD database and the 
VizieR catalog access tool, CDS, Strasbourg, France. S.T. is financially supported by the Japan 
Society for the Promotion of Science. 
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